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Abstract. We report the recent results of a search for lepton- flavor-violating r decays and a 
search for CP violation in r to uK^tv using a large data sample accumulated with the Belle 
detector at the KEKB asymmetric-energy e'^e' collider. The sensitivity to these modes is 
significantly improved compared to previous experiments. 



1. Introduction 

The success of the B-factory provides us the world-largest data sample of r+r^-pairs, that is 
useful for a search for New Physics (NP). Here, we discuss new phenomena relating to a r decays, 
in particular, lepton- flavor-violating r decays and CP violation in a r decay using an around 1 
ab~^ data sample accumulated with the Belle detector at the KEKB asymmetric-energy e'^e~ 
collider, 

An observation of the charged-lepton-flavor violation (cLFV) is a clear signature of NP since 
the cLFV has a negligibly small probability in the Standard Model (SM) even taking into account 
neutrino oscillations and many of models for NP naturally predict cLFV. In particular, lepton- 
flavor-violating (LFV) r decays are allowed to have various kinds of final states and some of 
them are expected to have a measurably large branching fraction, depending on the models. 
Therefore, not only to find the NP phenomenon but to distinguish the models, it is important 
to search for the LFV r decays. 

Similarly to the LFV r decays, an observation of CP violation (CPV) in a r decay is also 
a clear signature of NP. Here, we focus on the CPV in a r decay into uKgir induced by a new 
scalar particle such as a charged Higgs in a multi-Higgs-doublet model. This kind of CPV can be 
observed as a difference in the decay angular distributions by the interference effect between 
the scalar and (axial-)vector components of the intermediate state in the decay. 

2. T ^ £M° (i = e,/i, M° = 7r°,T7,T7',p°,K*°,K*°,a; and cf)) 

Generally, r fij is the most promising among various LFV r decays in many models for NP. 
However, in some models, such as non-universal-Higgs-mass model [1], r — )■ firj can have the 
largest branching fraction. In addition, a magnitude relation among the branching fractions for 
the LFV r decays will help to distinguish the NP models. Here, we show the results for the 
two-body r decays into a lepton and a neutral meson (pseudoscalar and vector). 



2.1. Analysis Method 

In the LFV r analysis, in order to evaluate the number of signal events, two independent variables 
are defined, that are signal-reconstructed mass and energy in the center-of-mass (CM) frame 



from energies and momenta for the signal r daughters. In the t ^r] case, they are defined as 



M^^ = JEI„-P^r„ (1) 



A 



where E^^^ {P^iri) is a sum of the energies (a magnitude of a vector sum of the momenta) for 
/U and r/, the superscript CM indicates that the variable is defined in the CM frame and the 
-^beam ™eans the initial beam energy in the CM frame. Principally, M^^ and /S.E are expected 
to be nir (~ 1.78 GeV/c^) and (GeV), respectively, for signal events while M^^ and AE' will 
smoothly vary without any special structure in the background (BG) events. The signal MC 
distribution of t ^ firj on the M^^ - AE plane is shown in Fig. [TJ Due to the resolution, the 




Figure 1. Distribution for the r — )• ;ur/ signal events. An elliptical region is defined as a signal 
region, that corresponds to ±3(7 and a region between two parallel lines excluding the signal 
region is used for the BG estimation as a sideband region, where the resolution a is evaluated 
by a asymmetric-Gaussian fit to the signal MC distribution. 

signal events are distributed around M^.^ ~ nir and AE ~ (GeV). Taking into account the 
resolution, we set the elliptic signal region. (See Fig. [TJ) Since the resolution for the photon 
energy is not so good as that for the momentum of the charged track, the signal region of the 
mode having some photons in the final state can be large. 

To avoid any bias for our analyses, we perform the blind analysis, i.e., we do not observe the 
data events in the signal region before fixing the selection criteria and the evaluation for the 
systematic uncertainties and the expected number of the BG events. Finally, we observe some 
signal events or evaluate an upper limit of the number of the signal event at 90% confidence 
level (CL) using the expected number of the BG event and the number of the observed events 
in the signal region by Feldman-Cousins approach. [2] 



Table 1. Summary for the efficiency (Eff.), the expected nmnber of the BG events, (N'^q) 
and the upper limit on the branching fraction (UL) for each mode, where "comb." means the 
combined result from subdecay modes. 
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2.2. T eP° (I = e, ^l, = ttO, r/, -q' ) 

We perform a new search for the r decay into a lepton (e or ^u) and a neutral pseudoscalar (vr*^, r/ 
or rj') with a 901 fb~^ data sample, that corresponds to 8.2 x 10^ r"'"r~-pairs. A neutral pion is 
reconstructed from 2 photons while an rj (?]') is reconstructed from 77 (p'^7) as well as Tr+vr^vr'^ 
(7r"'"7r~r7) to increase the detection efficiency. When the neutral pseudoscalars are reconstructed, 
their four-momentum is evaluated by a mass-constrained fit to obtain a better resolution for the 
signal region. Because we have modified the selection criteria applied to the previous analysis, 
we obtain an about 1.5 times better detection efficiency while a similar background level is kept. 
Here, the number of the BG events in the signal region is evaluated by an interpolation from 
the number of the observed events in the sideband region, assuming a flat distribution for the 
BG events in the signal and sideband region. 

As a result, we observe one event in the r — > e?7(— t- 77) mode while no events are found in other 
modes as shown in Fig. [2j Since these results are consistent with the background estimation, 
we set upper limits on the following branching fractions, as shown in Table [TJ Consequently, we 



(a)T-»ejt"(^Yy) 



O 



<l-0.4 ::.::r;::i::: 



0.2 




er|{r|^YY) 









0.2 














^ '^i. . 


0.4 




::;::»::::::: 








0.6 


: : ^'j 


• 


1.6 


1.7 1.8 



M,„o(GeV/tf) 



M (GeV/tf) 





0.2 









(Ge 


0.2 


m 


0.4 


< 






0.6 



(f)T-»|.l7t»(^ 


yy) 


. -/uU 








> . • •* 

■ i ^ , 



(g)T^liri(ri^Yy) 




. . . 











1.6 1.7 1.8 1.9 1.65 1.75 1.85 

A4,„(GeV/tf) M„„(GeV/tf) 



(c)T->er|(Ti-^Ti7t7i») : 



















0.1 - 
-: 
-0.1 ^ 
-0.2 -: 



• (d)T^eri'(ti'-»i'P°) 



(e)T->er|'(r|'-t.Ti:Ti:rii) 



-0.4 



1.7 1.75 1.8 1.85 

M^^(GeV/tf) 
(h)T ^ )iri ( r| ^ 71.7171:°;) 



4 1.78 1.82 1.72 1.76 1.8 

A4^.(GeV/tf) A4^.(GeV/tf) 



(i)T^Hti'(Tl'^Yp°) 



1.7 1.74 1.78 1.82 

A/, rGeV/c=) 



1.74 1.78 1.82 

A^,„,(GeV/tf) 



1.74 1.78 1.82 

A4,,,(GeV/tf) 



Figure 2. Resulting 2D plots for r — > e7r° (a), er\{-^ 77) (b), e?7(-± tttttt^) (c), e?7'(— > 7/5"(— > tttt)) (d), 
e77'(-> 7r7rr/(-^ 77)) (e), /xtt" (f), [L^{-^ 77) (g), [L^{-^ 7r7r7r°) (h), [irl {-^ 1P°{-^ tttt)) (i) and 
fj,ri'{^ TnTr]{^ 77)) (j), on the Mipo - AE plane. Here, black dots (shaded boxes) express the data 
(signal MC), the region bounded by two lines is defined as a 3a band for the BG estimation and the 
elliptic region is the signal one which corresponds to Scr in each plot. One event is found in the signal 
region for r — >■ ery while no events appear in any other modes. 



Table 2. Summary for the efficiency (Eff.), the expected number of the BG events, {N'^q) and 
the upper limit on the branching fraction (UL) for each mode. 
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obtain a (2.1 — 4.4) times more sensitive result than previously. 
2.3. T ^ eV° (£ = e,/x, F° = , K*^ , K*^ cjj) 

Similarly to — ?• i~P^, we update our results of the search for r" i~V^, where £ = e, /i, 
y° = p'^,K*^,K*°,uj,(l) with an 854 fb~^ data sample, that corresponds to 7.8 x 10^ r+r"- 
pairs |3j. Here, , K*^ , K*^ , co and (j) are reconstructed from tt'^'tt", tt~K^, tt^K^, Tr+vr^vr" 
and K~^K~, respectively. By performing a detailed background study, we obtain a 1.2 times 
better efficiency in average with keeping similar level backgrounds. For the modes including a 
fi, BG events mainly come from rr events while some two-photon process or radiative Bhabha 
process with a gamma conversion becomes main background in the modes having an electron. 
Finally, one event is found in the signal region for — t- fi^K*^, p,^ K*^ and p,^(j) while no 
events appear in any other modes. (See Fig. [3l) They are consistent with the expected number 
of the backgrounds. Consequently, we set the 90% confidence level upper limits on the branching 
fractions as shown in Table [2j At the present, the upper limit for r — )■ jxp^ is the most sensitive 
among the all LFV r decays. (See Fig. HI where all of the current upper limits are shown.) 
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Fi_gure 3. Resulting 2D plots for r ^ ep" (a), eK*° (b), eK*° (c), eu (d), e(f> (e), fip° (f), pK*° (g), 
pK*° (h), puj (i) and pcj) (j) on the M^yo - Ai^ plane. Here, black dots (shaded boxes) express the data 
(signal MC), the region bounded by two lines is defined as a 5a band for the BG estimation and the 
elliptic region is the signal one which corresponds to 3(j in each plot. One event is found in the signal 
region for t — >■ pK*^ , pK*^ and pcj) while no events appear in any other modes. 
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Figure 4. Upper limits on the branching fractions for the LFV r decays. Here, red, blue and 
black circles express the results obtained by Belle, BaBar and CLEO, respectively j^. 



3. CPV measurement with the r decay into uK^tt 

In SM, CP in the lepton sector is conserved. However, by introducing a new scalar field with a 
complex coupling constant, CPV can be induced. A charged Higgs in the multi-Higgs-doublet 
model [5J is one of the candidates for such a scalar field. Here, we discuss the CPV appearing 
in the r i^K^tt decay as an interference between the SM and NP processes [B]. 

3.1. Analysis method 

A hadronic current into K^n can be written as 

= {K^{qi)7r{q2)\u7^s\0) 

= {qi-q2r[g,u-^^^Fv{Q') + Q,FsiQ^), (3) 

where u and s are spinor fields for a u and s quark, respectively, qi and q2 are 4-momenta for 
and vr, respectively, Fy (Fs) is a vector (scalar) form factor and = Qi +Q2- To make the 
current induce CPV, instead of Fs, we introduce Fs, defined as 

Fs = Fs + ^Fh, (4) 
rrir 

where Fh is a form factor for the charged Higgs, r]s is a complex coupling constant [7]. The 
size of CPV is proportional to an imaginary part of tjs. From the current, in the hadronic rest 
frame {Q = 0), we can calculate the differential partial width for r vK'^tt as 

dr(T uK^ot:) 

— = (CP even terms) 



dio 



-4^\q-{\Q{FvF^)Q{r,s) cos V cos (3, (5) 

VQ 



where u = dQ'^d cos Od/S, 6 (/3) is an opening angle between the direction opposite to the one of 
the CM system and the one of the hadronic system in the r rest frame (between the direction of 
Kg and that for the CM system in the hadronic rest frame), and "0 denotes the angle between 
the direction of the CM frame and the direction of the r as seen from the hadronic rest frame. 
To extract the CP violating term, we define an asymmetry as 

,CP ^ /cos/3cosg(^-^)da; 

~ (cos/3cos^)^- — (cos/3cos'(/')r+- (6) 
Therefore, by measuring the second line, we can evaluate A'^^ experimentally. 

3.2. Result and Discussion 

With a 699 fb^^ data sample, we obtain around 3 x 10^ r vK^tt candidates, where it turns 
out that 23% of them are BG events by MC study. (See Fig. [^a).) Using them, A'^j^ is 

evaluated depending on \/Q^{= W). (See Table [3^ and Figs. El^b) and (c).) In the final result, 
'J — Z interference and vr^ detection asymmetry effects are corrected for and the background 
subtraction is also performed. 



Table 3. CP asymmetry in bins of the hadronic mass W. The first and second errors correspond 
to statistical and systematic ones, respectively. 





W (GeV/c^) 




0.625-0.890 


0.890-1.110 1.110-1.420 


1.420-1.775 




7.9 ±3.9 ±2.8 


1.8 ±2.1 ±1.4 -4.6 ±7.2 ±1.7 


-2.3 ± 19.1 ±5.5 



Consequently, no CPV is found. From this result, we can evaluate the 90% CL upper limit 
for |9(r7s')| varying a constant-strong-interaction-phase difference between Fy and Fg: 

|9=(r/s)[ < (0.012 - 0.026). (7) 

This result is almost one order of magnitude more restrictive than that previously obtained by 
CLEO [8J, i.e., \Q{t]s)\ < 0.19. 




Figure 5. (a) Distribution of W for the r — )• vK^tt candidate, (b) plot for depending on 
W and (c) zoom of (b). 



According to a theoretical prediction for Q{r]s) in the multi-Higgs-doublet model, rjs is written 

as 



where is an s-quark mass, Mjj± is the lightest charged Higgs mass, and X and Z are complex 
constants describing the coupling of the Higgs to the u and s quarks and the r and u, respectively. 
Using the limit < 0.026, we can obtain an inequality: 

|3(A-Z-)l<0.15^^. (9) 
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